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We consider a one-dimensional model for many-electron atoms 
in strong magnetic fields in which the Coulomb potential and in- 
teractions are replaced by one-dimensional regularizations associated 
with the lowest Landau level. For this model we show that the 
maximum number of electrons A'^max satisfies a bound of the form 
-^max < 2Z -|- 1 -|- c^/B where Z denotes the charge of the nucleus, 
B the field strength and c is a constant. We follow Lieb's strategy in 
which convexity plays a critical role. For the case N = 2 with frac- 
tional nuclear charge, we also discuss the critical value Zc at which 
the nuclear charge becomes too weak to bind two electrons. 
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1 Introduction 

It is well-known that systems in strong magnetic fields behave like systems in 
one-dimension, i.e., a strong magnetic field confines the particles to Landau 
orbits orthogonal to the field, leaving only their behavior in the direction of 
the field subject to significant influence by a static potential. Therefore, a 
better understanding of one-dimensional systems is essential to understand- 
ing the behavior of systems in strong magnetic fields. Although many one- 
dimensional systems, including that of a hydrogen atom with a single electron 
0, |], have been thoroughly studied, relatively little is known about multi- 
particle atoms confined to one dimension. 

In this paper we study the question of bounds on the maximum excess 
negative charge using one- dimensional models of many-electron atoms. Be- 
cause our goal is an understanding of the behavior in one-dimension, we do 
not deal with the question of accuracy of our one- dimensional models as ap- 
proximations to, or reductions from, real 3-dimensional atoms. However, we 
sketch such a reduction as motivation for the models considered. 

There is some question as to the proper replacement for the Coulomb 
potential in one dimension [§]. The potential V{x) = l/\x\ is so singular 
that the one-dimensional Hamiltonian —A — l/\x\ is not even essentially 
self-adjoint. Fortunately, an electron in a Landau orbit is attracted to a 
nucleus with charge Z, not by a potential of the form —Z/\x\, but by a reg- 
ularized potential which is finite at the nucleus. However, the corresponding 
regularization of the electron-electron interaction is more complicated unless 
both electrons have zero angular momentum in the direction of the field. In 
that case, the regularized interaction has a simple form which is the basis for 
our model. 

There is an extensive literature on atoms in magnetic fields. Interest in 
atoms in extremely strong fields, which began in the 1970's after the discov- 
ery of pulsars, has recently been renewed in the 1990's, culminating in the 
comprehensive work of Lieb, Solovej and Yngvason (LSY) |TB|, 0, 0, |2D| . 
For a discussion of early work on approximations for atoms in extremely 
strong magnetic fields we refer the reader to the insightful paper of Rau, 



Mueller, and Spruch [19| and to the review by Ruderman[^. References to 



later work are given in the introduction to |13| and a summary of the work 



of LSY |T^, |T8| is given in [|T^, ^. Rigorous work on atoms in magnetic 
fields began with the work of Avron, Herbst and Simon (AHS) [|l], |^, ^. 
LSY |TO, W7\, ITR B^] not only analyzed extensions of Thomas-Fermi the- 



ory in 5 distinct regions, but showed that these regions suffice to give the 
correct asymptotics for the exact Hamiltonian. In particular, they showed 
that the maximum number of electrons A^max(2^, -B) which can be bound to 
a nucleus of charge Z in a constant magnetic field of strength B satisfies 
lim ini NjT^s,x{Z,B)/Z > 2 as Z ^ cxd and B/Z^ —>■ oo. This should be com- 
pared with asymptotic neutrality [|15], |^ ^, i.e., lim.z-^ooN^i^^{Z,0)/Z = 1, 
for atoms without magnetic fields. We hope that the analysis of the simple 
model in this paper is a modest ffist step toward a better understanding of 
the mechanism by which extremely strong magnetic fields bind an "extra" Z 
electrons, as well as the conjectured converse iVmax(-Z', B) < 2Z. 

The full 3-dimensional Pauli Hamiltonian for an N-electron atom with 
nuclear charge Z in a constant magnetic field of strength B, acting on 7i^, 
the n-fold tensor product of L^(R ) (g) C^, is 
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H(iV, Z, B,a) = J2 [\Pj + Ap + a, . B - Z/\r,\] + J^ «/|r. " r.| (1) 

j=l j<k 

where A is a vector potential such that V x A = B and a is a cou- 
pling constant introduced for convenience in discussing scaling. If the spin- 
coupling term is omitted, it often suffices to consider the corresponding scalar 
Hamiltonian, which we denote H{N, Z, B, a), as an operator acting only on 
[L^(R^)]^ or the "space" portion of 7i^ . We will choose our coordinate 
system so that the field B = {B, 0, 0) is in the x-direction and the gauge 
so that 2 A = B X r. The Hamiltonian (1) satisfies the scaling relation 
li{N,Z,B,l) = BU{N,ZB-^/\l,B-^/^), i.e., we can scale out the field 
strength by replacing the nuclear charge Z by ZB~^^'^ and reducing the 
electron-electron interaction by B"^/"^. Alternatively, we could have included 
the electron charge unit e explicitly and replaced e by eS"^'^. 

We let Eq{N, Z, B, a) denote the infimum of the spectrum of the scalar 
Hamiltonian H{N, Z, B, a) defined above. It is well-known that the spectrum 
of the "free" Hamiltonian \Pj + Ap is [-B,oo], that the spectrum of |Pj + 
Ap + cTj ■ B = [cxj ■ (Pj + A)]^ is [0, oo) and that the lowest Landau level has 
energy B with infinite degeneracy indexed hj m = 0, 1, 2, 3, .... corresponding 
to angular momentum —m quantized in the field direction. Therefore, the 
continuous spectrum oi H{N, Z, B, a) is [B + Eo{N — 1, Z, B, a), cxd), i.e., the 
continuum begins at B plus the ground state energy for A^— 1 electrons. Thus 
the question of whether or not H{N, Z, B, a) has a bound state is determined 



by whether or not some test function \1/ satisfies 

(^, H{N, Z, B, a)^) <[B + Eo{N - 1, Z, B, a)] ||^||l (2) 

For the full Hamiltonian H(A^, Z, B, a), the continuous spectrum is [Eq^N — 
1,Z,B, a),oo) which implies that ii{N,Z,B,a) has a bound state if and 
only if there is a 'i' G 7i^ for which 

(*,H(A^,Z,5,a)*) <Eo(iV-l,Z,5,a)||*||^ (3) 

For the remainder of this paper we will omit explicit consideration of spin 
[although we will be able to draw some conclusions indirectly. See the remark 
after equation (|T3|).] Our methods cannot handle the a ■ B term explicitly 
and the inclusion of spin in the wave function does not affect the remaining 
results in any essential way. 

We now consider the the 3-dimensional energy minimization problem re- 
stricted to functions whose behavior orthogonal to the field is described en- 
tirely by product functions in which all electrons are confined to the lowest 
Landau level, i.e we restrict to N-electron functions \1/ of the form \l/r 



- mi...mjv 



<l>(xi ...Xn) Uk=i 7mfe in, Ok) where 



TTlfe' 



[7rm!]"^/^fi("'+^)/V™e"'^''''/^e"''^^ 



denotes the Landau level with energy B and angular momentum —m. (Note 
that we are using cylindrical coordinates {x, r, 9) with r = y^?/^ + z^ so that 
|r| = y/x"^ + r^.) Such expectations satisfy 

{^mr...mr,,H{N, Z, B, 1)^,„,...„^) 

= {^^„„^^,H{N,Z,B,l)<i/,n,...m^) + NB (4) 

= S(*™,...„^, HiN, ZB-^'\ 1, S-i/2)^^^.,.„^) + NB 

where 

N 
j=l j<k 

For each fixed choice of mi . . . rriN- the minimization problem can be reduced 
to a one- dimensional problem for the Hamiltonian H!^^''''^'^{N, Z,B,a) in 
which the kinetic energy has the usual — ^ form and the potentials l/|rj| 
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and l/\vj — Yk\ are replaced by regularized potentials denoted by V^. {xj) and 



^i^,n^S\^j~^k\) respectively. 



Vi{x) 



\-im?/\Y\dydz 

jDm+l poo 2mp — Br^ 

"ml Jo \/x^ + r2 



rdr 



\m\ 



11-1 



25 



m+l 



^Bx^ 



vn\ 



x^ 

oo 



ujB 



-.du 



(6) 



2\m^-BV 



it' - x'Te 



dt 



\x\ 



which satisfies the scaling relation V^{x) = \fB V^(\/5a;) as one would 
expect from the scaling properties of -ff(A^, Z^ B, a). W^.,^^ is defined anal- 
ogously; we postpone discussion of its explicit form. That the regularizations 
Vm,B{x) are important for atoms in magnetic fields goes back at least to Schiff 



and Snyder pj] in (1939) and played an important role in the AHS study 
of hydrogen. In the case m = and -8 = 1, (6) can be rewritten as 



l/o(x) 



v^e"'[l-erf(x)] 



(7) 



-.du = 2e' 



'dt 



'o Vx^ + u 

A comparison with the usual Coulomb potential is given in Figure |I]. 

We now describe the correspondence between a restricted three-dimensional 
minimization problem for a Hamiltonian with a magnetic field, namely. 



7-imi...mjv 
^0 



{N,Z,B,a)= inf (^^,...^^, i7(iV, Z, S, «)^^,..^^) (8) 



* 



and the one-dimensional minimization problem 



77imi...mjv 
^0,x 



{N, Z, B, a) = inf ($, H^^-'^n ^p^^ z, B, a)*) 



where 



N 



J=l 



d"^ 



j<k 



mk\\^j 



(9) 

Xk\) 

(10) 




Figure 1: Comparison of the potentials (reading from the top down) — ,} , 
-Vo{x), and -i^i 

in which no magnetic field is explicitly present. Since 

{^m,...n.j,,H{N,Z,B,a)^m,...n.^) = {<!> , H:''-"''' {N, Z, B , a)<!>) + N B (11) 
we can conclude that 
(^ 



- mi...mjv/ ^ ^0 



mi...mjv-i 



:^,,„^^,H{N,Z,B,a)^ 

^^ (<l>,ifr-™-(iV,Z,5,a)<l>) > Eo" 



(A^- l,Z,B,a) + B 
"''■■■"'''-\N-l,Z,B,a). (12) 



Thus, we have reduced the problem of determining whether or not the some- 
what artificial problem of whether or not an atom in magnetic field B whose 
electrons have prescribed angular momentum corresponding to mi . . . tjin has 
a bound state to that of whether or not a related one-dimensional system 
has a bound state. Although, we do not consider spin-coupling explicitly, it 
does follow from (3) that we can also conclude that 



(*„,...„^,H(Ar,Z,5,a)* 
^^ ($,^™i-"'^(iV,Z,5,a)$) > El 



m,...m^)>El'-"'--^{N-l,Z,B,a) 



mi...m]v-i 

X 



(N 



l,Z,B,a). (13) 



where '^mi...mN is chosen so that all components in C^ correspond to spin 
"down" . (This observation, which is a fortuitous consequence of the physical 
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fact that the coefficient of a ■ B is exactly 1, was pointed out to us by J. P. 
Solovej.) 

We now consider the regularization of the interaction in the special case in 
which mj = m^ = 0. Recall that 7^(1/, z) = Bn-^e-^^y^+'^^l^ = Bir-^e-^''^/'^. 

rBn^ . ,^ [[,..,. f [,..,. ITo'^IZ/I, ^l) I' l7<f (2/2, ^2) P 



Wq^q{\xi - X2\) = dyidzi / / dy2dz2 



a/(xi -X2)2 + \{yi,zi) - iy2,Z2W 
B / rlf ^ (14) 

Jo V(Xi - X2)2 + 2t2 

_]_yB f\Xl-X2\ 



where we have made the change of variables 

s = ^[(2/1' ^1) + (2/2. ^2)], t = -7^[{yi, Zi) - {y2, Z2)] 
with s = |s|, t = |t|, and used the fact that 

b^iyuz,)\'b^iy2,z2)\' = 5V-V^(^'?-^^i) 

= 5^7r-V^(^^+*^) = |7o^(s)n7o^(t)p. 

Thus, the exceedingly simple relation iyQ^Q(|a;i— X2I) = 2^^^'^Vf^ {2~^^'^\xi — X2I) 
follows from the invariance of |7(f (s)p|7(f (t)p under the transformation of 
s,t to (s ± t)/A/2. This unusual invariance, corresponding to the mixing 
of coordinates of two particles, will not hold if m 7^ 0. Symmetrizing the 
product or replacing \l/mi...mjv by an arbitrary element of the projection onto 
the lowest Landau level, would require consideration of exchange terms as 
well. Therefore, we will only study models corresponding to constraining 
all electrons to have angular momentum zero. Since \l/o...o is then symmet- 
ric with respect to exchange of {yj,Zj) with {yk,Zk), it will have the same 
permutational symmetry as $(xi . . . xn)- 

Despite the severity of the restriction to m = 0, our model seems well- 
suited to study of the bounds on the negative ionization. Integration by parts 
of (^ easily yields 

V^^.ix) < V^ix) < V.^ix) < l/\x\ (15) 

with the difference greatest at the origin and all potentials satisfying V^{x) ~ 
l/\x\ for large x. Moreover, for any choice of rrij, rrik whenever \xj — Xj\ is 



large, Wj^.^^^{\xj — Xj\) ~ l/l^j — Xj\ as well. Thus it appears unlikely that 
placing some electrons in Landau levels with m ^ will allow binding if 
m = does not. 

With this heuristic background, we study models for N-electron atoms in 
one-dimension corresponding to Hamiltonians of the form 

h{N,Z,M) = MH^-\N,Z,M-^,1) (16) 

N 



i:[-fI^-^H^5X^) 



j<k 

with the "mass" M proportional to B^^/"^. Because of the scaling relation 
H{N, Z, B, 1) = BH{N, ZB~^/^, 1, B~^/^), the only role of the field strength 
in the one- dimensional situation, is to reduce the mass by a factor of B^^'"^. 
Thus, for simplicity, we have set both B = 1 and the coupling constant a = 1 
leaving the field strength implicit in the mass M. 

As observed in 0, the regularized potentials satisfy 

^ < ^ = < V^{x) < ^. (17) 

^{m + l)B + \x\ ^{m+l)B+ |x|2 \A 

The cut-off potential Kut(a^) = l/(kl + 1) is also of some interest. Haines 
and Roberts |^ have given explicit solutions to the eigenvalue problem for 
the hydrogenic Hamiltonian —A — Kut(a;). By the above remark, V^ut(a;) < 
V,\x) < l/\x\. 

In sections 2 and 3, we study Hamiltonians of the form (0) in the two 
cases V = V^ut and V = Vq. For both models we show that the existence of 
a bound state implies A^ < 2Z -|- 1 -|- c^/B. 



2 Two-electron systems 

We now discuss in more detail the behavior of the one-dimensional Hamilto- 
nian ( p!6D when the number of electrons is iV = 2. Although our discussion 
is descriptive and non-rigorous, we believe the insights are useful. The two- 
electron Hamiltonian can be written in the form. 



/i(2,Z, 5-1/2) = -V5 






+ W{xi,X2) (18) 



where 

W{x,, X,) = -ZV,{x^) - ZVo{x2) + -^Vo (^^^^) (19) 

Consider a classical system of three particles on a line with charges —1, 
+Z, —1 interacting with the usual inverse square Coulomb force (or, equiv- 
alently, qiqj/\xi — Xj\ potential). The only potentially stable configuration 
is one in which the particles are arranged symmetrically with the two nega- 
tively charged particles on opposite sides of, and the same distance from, the 
positive "nucleus" as shown in Fig. 1. Even in this case, the system is stable 
only for Z = 1/4. For Z < 1/4, both "electrons" will move off to infinity, 
while for Z > 1/4, the "electrons" collapse into the center. This suggests 
that for Z > 1/4 one might be able to show that the system binds by using a 
trial function "^a in which the two electrons are localized on opposite sides of 
the nucleus at a distance a sufficiently far from the center that Vo(±a) ^ 1/a. 
(Permutational symmetry is not relevant; if such a trial function binds, then 
a bosonic system will also bind.) This picture would require that Axi < a 
for each electron so that, by the uncertainty principle, the kinetic energy 
satisfies {"^ a[- -£t]^ a) > l/a'^- Thus we estimate 

(vl/,/.(2, Z, B-'/')^,) > 2^ - — + ^. (20) 

a^ a 2a 



Minimizing over a yields amin = 2\fB jiZ — i) and 

(21) 



{^ji(2,z,B-^i^)^:)>-^^ '^' 



2/B 

The continuous spectrum begins at the ground state energy of the corre- 
sponding one-electron Hamiltonian /i(l, Z, B^^/'^) = —yB^ — ZVo{x) which 
Avron, Herbst and Simon p| showed is given asymptotically (for large B) by 

Eq{1, Z, B^^/'^) = —4^ (log 4^ j Thus, binding would require 

{Z-\r>2Z^Lg^y. (22) 



This is obviously false for large B, which suggests that the uncertainty prin- 
ciple prevents the corresponding one-dimensional system from binding, even 





Figure 2: W{x, —x) and W{x, 0) for (reading from the top down) 
Z = 0.25, 0.3, 0.35, 0.4, 0.5, 0.6, 0.7, 0.8. The left graph, W{x, -x), describes 
the potential when the electrons are on opposite sides of the nucleus with 
each a distance x away. The right graph, W{x,0), describes the potential 
when one electron is fixed at the origin; the left asymptote then corresponds 
to the other electron moving to infinity. 

when the classical system collapses. Because M ~ 5~^'^, large field strength 
corresponds to small mass, i.e., as the field strength increases, the electrons 
become "lighter" and more difficult to localize. Thus we must seek a different 
mechanism to explain binding of a one-dimensional two-electron atom. 

Fortunately, our model provides an alternative mechanism for binding. 
The regularization of the potential at the origin combined with the effective 
reduction of the interaction by l/v2, permits both electrons to be close to 
the nucleus for Z > 1/2^/2 in the sense that -2ZVo{0) + T^I'^Vq < 0. Of 
course, the uncertainty principle also precludes binding with a trial function 
in which both electrons are exactly at the center. Nevertheless, we believe 
that the mechanism for binding is that the effective reduction in the interac- 
tion permits the two electrons to overlap strongly near the nucleus. 

The Hamiltonian (0) can be regarded as describing either two particles 
in one-dimension or one particle in the field of the two-dimensional potential 
W{x,y) given by (0). Regardless of our viewpoint, the continuous spec- 
trum will begin at the AHS estimate of —4^ x^'^^'Tr) ■ Since this is 
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close to zero for large B, we will regard W{x, y) as attractive where it is 
negative and repulsive where it is positive. One might expect W{x,y) to 
have its minimum on the line y = —x, corresponding to electrons on op- 
posite sides of the nucleus. The actual situation, which is described below, 
is more complex. Plots of W{x, —x) and W{x, 0) are shown in Fig. 2 for 
Z = 0.25, 0.3, 0.35 ^ ^, 0.4, 0.5, 0.6, 0.7 ^ ^, 0.8. The graph of W{x, 0) is 
deliberately asymmetric so that the left asymptote, W^(— oo, 0) can be com- 
pared to the minimum of W{x, —x) at a; = 0. We now discuss this behavior 
in more detail, noting that all the qualitative features can also be verified 
analytically. 

• J < Z < -j= : The potential does have a pair of weak minima along the 
line y = — x; however, closer examination of the full two-dimensional 
potential shows that these are not true minima, but saddle points for 
W{x,y). The potential is repulsive at the origin and only weakly at- 
tractive elsewhere. 

• -K= < Z < ^ : The potential is now attractive at the origin. However, 
as above, the two weak minima on the line y = —x correspond to saddle 
points of W{x,y). 

• \ < Z < -j= : The two weak minima on the line y = — x have now 
coalesced into a true (two-dimensional) minimum at the origin. How- 
ever, this is only a local minimum. W{0, 0) > W{0, oo) so that the 
energy will decrease if one of the particles remains at the origin while 
the other goes off to infinity. Fig. 2 suggests that the minimum along 
the X = and y = lines is too shallow to prevent one of the electrons 
from tunnelling through to infinity. Thus, we expect resonances, but 
not bound states in this region. 

• Z > 4= : The potential has a true minimum at the origin. 

Thus, the behavior of the potential W{x,y) strongly suggests that, at 
least for sufficiently large B, binding occurs for Z > 4= Ri 0.7 and that the 
mechanism which permits this is the reduction in strength of the repulsion 
by l/v2 which permits both electrons to simultaneously sit near the nucleus. 
The behavior of the regularized potential model for a two-electron atom in 
one dimension seems to be very different from that of classical Coulomb 
particles confined to a line. The binding cut-off of Zc ^ 0.7 represents a very 
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Figure 3: Saddle point of W{x, y) for Z = 0.4 in the quadrant x < 0, y > 0. 

modest "binding enhancement" when compared to the value of Zc ~ 0.9 
obtained in |^ for a two-electron atom in three dimension in the absence of 
a magnetic field. [An "enhancement" to the level of 2 = A^ > 2Z + 1 would 
require Zc < O.5.] 



3 Kinetic Energy and Lieb's Strategy 



We now describe the elegant strategy used by Lieb |jT3| (see also section 3.8 
of [^]) to show that for real atoms the maximum number of electrons that can 
be bound is less than 2Z + 1 irrespective of permutational symmetry. The 
essence of Lieb's argument is to show that one can replace the variation over 
N-electron wave functions \E'(ri, . . . r^r) by a variation over density matrices 
of the form 



r(ri, . . . Tjv; Si, . . . S]\r) 



(23) 



N 






where z/(r) > is a strictly positive function which will be chosen later. 

1/2 

Note that, even though the functions \l'z/- are not orthonormal, TrFA = 
^ .([z/(rj)]^/^\E', A[z/(rj)]^/^\l') for any operator A. Despite the introduction 
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of the function v, which will be chosen later, the \& used in F satisfies the 
same symmetry or domain constraints as the original minimization problem. 
It will be useful to describe Lieb's argument in the rather general situation 
of an N-particle Hamiltonian, Hj^ with the structure H^ = H-'^_^ + Kj + V^^^ 
where Kj = \Pj + Ap denotes the kinetic energy, H-j^_^ is the A^ — 1-electron 
Hamiltonian in which the j-th electron is omitted and V:;^^ is a potential which 
describes the interaction of the j-th electron with the rest of the system. Lieb 
studied H^ = H(N, Z, 0, 1) while we will consider Hn = h{N, Z, M). In our 
case, 

VL = -ZV{x,) + J]2-VV(2-V2|^^. _ ^^1) (24) 

Now assume that if at actually has a bound state so that there is a \&o 
satisfying Hn^Iq = -Eq^o- This implies that for any function u and any choice 
of J, 

z/(r,)i7jv^o(ri, • • • r^) = Eoiy{rj)'^o{ri, ...r^) (25) 

as well. Hence 

(^0, i^{r,)HN'^o) = Eoi^o, ^^(r,)^o )• (26) 

Then if Fq is the density matrix of the above form corresponding to the 
ground state \l/o, and Uj denotes i^{rj) 

N 



Eo{N)TiTo = Eo{N) J](^z/,^) 



N 

> EoiN - l)TrFo + J^ [i^f"^^ ^ ^T"^) + (^'^.^.^) 



The inequality in ( [27| ) above follows from the variational principle for ii\^_^. 
Although I^VjKj'^) need not be real in general, a careful analysis of the 
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argument above shows that when \1/ is an eigenstate of if at it is real because 
it can be written as the difference of two real quantities. In Lieb's original 
formulation, v was chosen so that this real quantity satisfied {^VjKj^) > 0. 

Then if Ej>J^^^^>j^^^) > it follows that Eo{N)TtTo > Eo{N-1)TtTo 
which contradicts the requirement for binding of Eq{N) < Eq{N — 1). This 
reduces the problem of showing that H{N, Z, 0, a) has no bound states in the 



1/2, 



1/2, 



absence of magnetic fields to that of showing that ^Aj^j "^V^f-Uj' \E') > 0. 
However, in a constant magnetic field the continuum begins at B + Eq{N — 1) 
so that a similar reduction would require the stronger condition {^VjKj^>) > 
B. 

In order to study (\l'z/jiC,\E') in more detail in different situations, it is 
convenient to write z/j as g"^ and observe that formally 



(V, \tVj + Apvl/) = ([^V, + A]{g^)g, [zV, + A]{g^il)g-') 



([zV, + A]g^, [zV, + A]ig^)) - J \g^\' 
±2i^{^Vg,[iVj+A]{g^).) 



^9 



(28) 



Note that the localization error, — J |(7\E'|" 
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/|(7^pV(7-V((7-i), shows 



a certain symmetry between g and g ^. If we now let g 



V 



1/2 



^{^VjK^^) = (^z/, 



1/2 



K,- 



^0^0 



2v, 



^-T) 



(p8|) implies 



(29) 



For the typical choice v = [V{r)] ^, one has Vz^(r)/z/(r) = — W(r)/y(r) 
which yields 



^{'^UjKj'^) = {^u. 



1/2 



K, 



V,V, 



2V, 



^uy^) 



(30) 



in accordance with the symmetry between g and g^^ noted above. Because, 
as discussed above, we will only be concerned with applications for which 
(A^VjKj^)) is real we will henceforth omit the 3fJ. We now discuss several 
cases in more detail using d to denote the dimension of the space on which 
Kj acts. (The original proof of Lieb used an argument which originated 
with Benguria (see Lemma 7.20 of [|l^) in the spherically symmetric case 
to show directly that 9fJ(z/\E'(— A)\l') > 0. The variant given here is due to 
Hoffman-Ostenhof 0. For other strategies see [^, ||, |10| ) 
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• A = 0, d > 3 : In this case we are interested in z/ = [V^(r)] ^ for poten- 
tials, (particularly the usual Coulomb potential V^(r) = l/|r|) which satisfy 



vy(r) 

V{r) 



< ^. Then K 



-A and (|30| ) becomes 



> {^i^. 



1/2 



-A, 



4|r| 



^z/'/') > 



(31) 



» A 7^ 0, d = 3 : In this case we can only conclude from (pODthat 

1/2 II 2 /,T, 1/2 



{^UjKj"^) > 5||^z/, 



{^-r 



dx] 


2u, 


2 



1/2 



^v]^') (32) 



> BW^v 



1/2||2 



4z/,- ' 



We could instead have proceeded as in ( pi]) above to obtain, (\E'i/j/('j\E') > 0, 
but we cannot conclude that i^VjKj^) > i?||\E'z/- p. The kinetic energy 
in the field direction is not able to control a 3-dimensional potential, 



Pi 



92 



such as 1/r^ arising from 



2v, 



Our one- dimensional models circumvents 



this problem because the choice of product state involving Landau functions 
ensured that the Py + Pf terms took care of the B, leaving P^ to deal with 
a one- dimensional potential[\^'(a:)/2V^(a;)]^ as described below. 

• A = 0, d=l: In one dimension, A plays no role. Although we will still 

2 



< 



2 it is not 



choose u = [V{\r\)] ^ with potentials which satisfy 

true in one or two- dimensions that —A — |2x|~^ > 0. Instead we will treat 



V'jx) 



\x\ 



\YL\ 

\2V\ 



as a potential in 



(v^z/,(-A,)vI/) = (M/z/; 



1/2 



d^ 


V'{x,) 
2V{x,) 


2' 



^uf). 



(33) 



Now in one-dimension the everywhere negative potential — |||7 1 will always 

Thus we can only conclude 



nve rise to a bound state of — -rr — 



that (^z/j(-Aj)^) > -e||^, 
operator. However, the bound 



1/2||2 



dx^ 



2V\ 



where — e is the lowest eigenvalue of this 



(^z/j(-Aj)^) > -cj||^| 



(34) 
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where uj = sup^. |z/'(x)p/4z/(x) will be more useful. Notice that, unlike the 
3- dimensional case where the lower bound is zero, when applying this result 
to if AT = h{N, Z, M) we will need to take into account the fact the entire 
kinetic energy term is multiplied by 1/M (or y/B). 

4 Bound on the Maximum Negative Ioniza- 
tion 

We now apply Lieb's strategy, which yields N^^a^iZ, 0) < 2Z + 1 for atoms 
without magnetic fields, to our one- dimensional models for systems in strong 
magnetic fields. We obtain the following result. 

Theorem 1 The maximum num,her of electrons A^max for which a Hamil- 
tonian h{N, Z, B^^^"^) of the form ^Tdj ) has a bound state satisfies A^max < 
2Z + 1 + c\^ for some constant c. 



In the interesting case B = 0{Z^) (which is the boundary of the LSY hyper- 
strong limit region in |T^), this yields a bound of the form A^max < 2Z+cZ^^'^, 
rather than a linear one of the form A^max < cZ or the expected optimal 
Ar^ax<2Z + o(Z). 

We apply the strategy of Section 3 with H^^ = h{N, Z, M) and u = 1/V. 
In this case, the analysis of expectations of V^^^ is straightforward because 

= -Z ^ ||M/f + ^ 5^(M/, 2-VV(2-V2|x, - x,\)V{x,r'^) (35) 

j j k^j 

= -NZW^W' + 5^(vl/, 2-i/V(2-i/2|^, - x,\)[Vix,)-' + Vix,)-']^). 
j<k 

Thus, if the potential satisfies 

2-VV(2-V2|a;^. _ x,\)[V{x,)-' + V{x,)-'] > 1, (36) 

then 

E^'^f *^nt^f ^) > [-NZ + N{N - l)/2]m'. (37) 
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For the cut-off potential V{x) = l/{\x\ + 1), ( PU[ ) follows easily from the 
triangle inequality since 

2'/'[V{2-'/^\w-x\)]-' = \w-x\+2'/^ 

< \w\ + 2^/2 + |a;| + 2^/2 = [V{w)]-^ + [V{x)-'] 

For the regularized potential Vo{x) we will use instead the convexity of 
i'o{x) = [Vo{x)]~^. (Because the proof |^ ^ of this essential fact, although 
elementary, is rather delicate and not readily accessible, we provide a sketch 
in the Appendix.) Thus we find, 

[Vo{w)]-' + [Voix)-'] > 2[Voi\w - x\/2)]-' 

> 2^/^[Vo{2-^/^\w - x\]-^ (38) 

where the first inequality used Vo{x) = Vo(|x|) as well as the convexity of 
z/o(a;) = [Vo(x)]^^ and the second inequality follows from 

Vm{2-^l^y) = [ml]-' / ^=^ dt 

Jo 






fm -t 



< 2'/'Vmiy) 



with y = 2-'^'^\w — x\ and m = 0. 

Thus we can conclude that, for both the cut-off potential VJ;ut(a;) = 
l/dsl + 1) and the regularized potential Vo{x), the one dimensional Hamil- 
tonian h{N, Z, M) satisfies 

Trro/i(Ar, Z, M) (39) 

> Eo{N, Z, M)Trro + [-NZ + N{N - l)/2 - Nuj/M]\\m\\^. 

The second term in ( ppj ) will be positive if (A^ — l)/2 > Z + uo/M. For both 
potentials, u = sup^. |i/'(x)p/4z/(x) is given by liinx^o\i''{x)\'^/4:h'{x). This 
yields, Wcut = 0.25 for the cut-off potential and coq = vr"^/^ < 0.18 for the 
regularized potential Vo{x). The latter follows from uj = z/(0)'^ = 7r~^/^ and 
the fact that, as shown in the Appendix, |z/'(a;)p/4z/(a:;) = z/(x)[z/(a;) — |x|]^ 
is decreasing for a; > 0. Thus, in both cases, the second term in ( pOf ) will 
be positive if A^ > 2Z + 1 + 1/2M. Since binding implies TrToh{N, Z, M) < 
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Eq(N, Z , M)TtTq and M is proportional to B ^/^, we have shown that one 
can find a constant c such that our model one-dimensional system does not 
have bound states if iV > 2Z + 1 + cv^ or A^^ax < 2Z + 1 + cVB. 

Although this bound is not optimal, it should be remembered that we are 
analyzing a Hamiltonian with the structure 

h{N, Z, 5-1/2) = h{N-l,Z, B-^'^) + B^/^K, + V^,{xj) (40) 



and that the factor y B multiplying the kinetic energy is something of a 
two-edged sword. On the one hand, it raises the energy of the effective 
one-electron Hamiltonian B^/'^Ki + J Vint\i^{xi^X2, ■ ■ ■ , XN)\'^dx2 ■ ■ ■ dxN] on 
the other, it multiplies any error arising from the kinetic energy - whether 



using Lieb's strategy or the Ruskai-Sigal localization approach 0, |21|, ^ |2^ 
- by v5- Since such correction terms are typically negative, the result of 
such treatments is to perversely magnify the negative error associated with 
a positive term otherwise expected to raise the energy. 
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A Convexity and Properties of l/Vo(x) 

Because ^{x) = l/Vo{x) is nearly linear (see Figure ^), the proof of its 
convexity is somewhat delicate. It will follow from the upper bound in the 
following pair of inequalities, which are of some interest in their own right. 



< u[x) < . (41) 



4 ' ' 3 

We now restrict attention to x > and define 



k 

9k{x) = , (42) 

^ ^ {k-l)x + V^^Tk ^ ' 
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Figure 4: Comparison of l/Vo(x) and x for x > 0. 



so that (^TD is equivalent to 

Qiix) > Vo{x) > g^ix) 



(43) 



which is an improvement on the classical inequalities of Komatsu |11|. The 



upper bound and convexity were established independently by Wirth ||2^ 



and by Szarek and Werner 1^^; the proof of lower bound which follows was 
communicated to the authors by E. Werner. 

We now observe that Vo(0) = y/n and Vq satisfies the differential equation 

(44) 



V^{x) = 2[xyo(x) - 1] 

and look for analogous behavior for g^. Since g[{x) = 
and xgk{x) — 1 = - — t\ gk{x), one can verify that 

g'k{x) > 2[xgk{x) - 1] 



k {k-l)Vx^ + k + 



X 



< 



2k 



y/x'^ + k {k — l)x + \/x^ + k X + \/x'^ + k 
(k - 2)x2 + k{k -3) <{k- 2)xVx^ + k 
x\k-2){k~4) + k{k-3Y <0. 



(45) 



For A; = 3, we use the first equivalence to conclude from x^ < xy/x^ + 3 that 
g'^ix) > 2[xg3{x) — 1]; for /c = 4, we use the second to conclude k{k — 3)^ > 
=> g'^{x) < 2[xg4ix) - 1]. 



19 




Figure 5: Graphs showing (74 (x) > Vo(a;) > (73 (x) for x > 0. 



For fc = 3, let /i3(x) = Vo(x) — gz{x). It is easily to verify that /i3(0) = 
Tx — -\/3 > and that \m\x~*oo \h2,{x)\ = 0. The first equivalence in (^Sf) 
implies that x^ < xvxM-3 =^ g'si^) > 2[xf73(x) — 1] so that h'^{x) < 
2x[Vo(x) — f73(x)] = 2x/i3(x) Vx. Now suppose that for some x = a, /13(a) < 0. 
Then h'^{a) < and it follows that /i3(x) is decreasing for x > a, which 
contradicts limj.^oo ^3(2;) = 0. This proves the lower bound in (^). The 
upper bound is proved similary ||2^ except that one now shows that h^i^x) = 
gi{x) — Vo(x) > 0. It is interesting to note that the upper bound in ( p3D is 
optimal, but the lower bound can be improved [^ to g-n{x) < Vo(x). The 
tightness of these bounds is evident in Figure |^. 

In order to use this to establish the convexity of v{x), we note (|4^) implies 
that p{x) = [Vo(x)]~^ satisfies 



Then 



z/'(x) = 2z/(x)[z/(x) — x]. 



iy"{x) = Aw' - 2u'x -1v = 4z/(z/ - x)(2z/ - x) 



2v 



ll^(^' + 4)] 



%v\(v - |x)' 

8z/ {{v{x) - lM(x) + 17^^+4)2 - (iv/^^T4)2] . 



(46) 



(47) 



Since ( ^iD implies z/(x) > \lg^{x\ we can conclude that v"{x) > 0. 

To show that |z/'(x)p/4i^(x) = i^(x)[i^(x) — |x|]2 is decreasing for x > 0, 
we note that (^Hf ) yields 

— [z/(z/ - x)^] = 2z/(z/ - x) [3z/2 - 4z/x + x^ - 1] 
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6z/(z/ — x) 



[V 






9 



<0. 



and (^rp implies that the expression in square brackets is negative. 

It is useful to have bounds on the extent to which v{x) and Vq{x) deviate 
from |x| and the Coulomb potential, respectively, when x is large. The upper 
bound in (p|) implies \i'{x) — \x\\ < gA{x) — x < l/2|x| while the lower bound 
yields |Vo(x) — l/|x|| < 1/x — gz{x) < \/2x^. The latter can also be proved 
directly from the integral representation (^ . 
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